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We present a calculation of 1/ml corrections to the lifetime differences of Bs mesons in the heavy- 
quark expansion. We find that they are small to significantly affect AF (Ba) and present the result 
for lifetime difference including non-perturbative 1/mb and I /ml corrections. We also analyze the 
generic AB — 1 New Physics contributions to the lifetime difference of Ba mesons and provide 
several examples 



Mixing phenomena in heavy bosons system is con- 
sidered as an important test of Standard Model(SM) 
and a probe for New Physics(NP) beyond it. Usu- 
ally it is referred to the fact that such process oc- 
curs only at the one loop level in SM . This makes 
it sensitive to the effects of new particles running in 
the loop. These interactions induce non-diagonal ele- 
ments in mass-matrix making flavor and mass eigen- 
states to be different. Analysis of mixing in charm, 
beauty systems led to positive signals which seem to 
be very well explained by Standard Model physics. 
The lifetime difference AP^ is generated by on-mass- 
shell intermediate states and seems to be one more 
test of Standard Model and heavy quark expansion. 
Yet some contribution of NP is still possible as an indi- 
rect probe of energy scales beyond currently accessible 
at experimental facilities. Further in paper we set up 
relevant formalism and discuss the need to compute 
corrections. After these corrections computed 
we consider impact of New Physics Ab — 1 interac- 
tion on the numerical value of AF. 



I. FORMALISM 

The width difference between mass eigenstates is 
then given by [l| 



Ar« =rr -Ff. = -2ri2 = -2F 



21, 
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where Fy are the elements of the decay-width matrix, 
= 1,2(11) = 15,), |2) = |B:)) 
Using optical theorem, off diagonal elements of mix- 
ing matrix can be related to the imaginary part of the 
forward scattering amplitude: 



2M 



(2) 



T = Im^y" d''xT{H^s{x)He«{0)}. (3) 

where ^^e// is an effective weak hamiltionian de- 
fined as follows: 



6 



(4) 



where four-quark operators are defined in the fol- 
lowing way: 

Ql = ihcj)v-A{cjSi)v-A, (5) 
Q2 = ihCi)v-AicjSj)v~A, (6) 

Q3 = {hsi)v~A{qjqj)v-A, (7) 
Q4 = (biSj)v-Aiqjqi)v-A, (8) 

Q5 = {hsi)v~A{q3q3)v+A, (9) 
Qe = (biSj)v-A{qjqi)v+A, (10) 
Qs = £^m,ka^'^{l-^5)Tt,s,G^,,. (11) 

In the heavy-quark limit the energy release is large 
and process is dominated by short-distance physics. 
An operator product expansion can be constructed 
which results in series of operators suppressed by pow- 
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ers of 1/mf,: 
T{Bs)2i 



2M 



l-J2{B4%\B, 



5:^(i?.|Or=^(^)|i3.). (12) 



The most recent calculations of Bg lifetime difference 
[l[ and of QCD corrections to AF^ do not provide 
definitive theoretical prediction of its value. Heavy 
quark expansion corrections of order of l/rub appear 
to be about 25% of leading order and QCD correc- 
tions are as big as 30%. [11]. We compute 1/m^ cor- 
rections in heavy quark expansion to directly check 
convergence of this series. In other words we com- 
pute matching coefficients of an effective Ab = 2 la- 
grangian. Computation of matrix elements of most of 
these operators is rather difhcult task due to lack of 
results from Lattice QCD and Light cone QCD calcu- 
lations. We used a factorization approach to estimate 
matrix elements of such operators. 

Expanding the operator product ([3|) for small x 
l/rUb, the transition operator T can be written, to 
leading order in the 1 / mf, expansion, as P, 0| 



Utt 

which results in Q 

r2i(B,) = -- 



iVclVj' [F{z)Q{^i2)+Fsiz)Qsit^2)] 



(13) 



-(v;;Ks)'VT^ X 



12tt{2MbJ 

X { [(1 - z) (2 C1C2 + N,Cl) + (1 - 4z)CiV2] (Q) 
+ {l + 2z){2C^C2+N,Cl~Cl){Qs)] , (14) 

where z — ml/ ml and the AB — 2 operators are as 
follows: 



Q = ibiSi)v-AibjSj)v-A, 
Qs = (biSi)s-pibj'^j)s-P ■ 



(15) 



Color re-arranged operators Q — {biSj)v-A{bjSi)v-A 
and Qs — {biSj)s-pibjSi)s-p that appear during cal- 
culations were eliminated using Fiertz identities and 
equation of motion. 

The Wilson coefficients F and Fs are obtained by 
computing the matrix elements of T in ([3]) between 
quark states. 

The coefficients in the transition operator at 
next-to-leading order, still neglecting the penguin sec- 
tor, can be written as [2]: 



Fiz) = Fniz)Cl{fii) + Fi2{z)Ci{fii)C2{tii) 



+ F22{z)Cl{ni), 



F,,{z)^f'^\z) 



as(A^i 
47r 



(16) 



(17) 



Fs{z) has similar structure. The leading order func- 
tions pjj^ , Fg'^j read explicitly 



Fi?{z 
F<"^ (z 

f(,")(z 
F^°^ (z 

F^^iz 

FfUz 



3Vl -4z(1 - z), 
3V1-4z(1 + 2z), 



= 2Vl -4z(1 - z) 
= 2\/l - 4z(1 + 2z), 

= i(l-4z)3/2 
2^ ' 

= -Vl-4z(l-H2z). 



(18) 
(19) 

(20) 



The next-to- leading order (NLO) QCD expressions of 
fI^\ F^gl- and corrections to Eq. (fT^ arrising from 
penguin diagram are given in Ref. [2] . 



II. 1/ml CORRECTIONS 

The general expression for lifetime difference of Bg 
mesons can be presented in the following way: 



F2i(S,) = — 



Ps{z)]Qs (21) 



127r(2A/Bj 
X {[F{z) + P{z)]Q+[Fs{z 

where 61 /m and 61 denote contribution from opera- 
tors suppressed as 1/mb and l/m^ respectively. These 
terms and their numerical values are computed fur- 
ther. 

The matrix elements for Q and Qs can be 
parametrized in the following way [l|, 0, 

{B,\Q\B,) = fl Ml2\l + ^]B, 



{Bs\Qs\Bs 



-fl Ml 



Ml 



(22) 



{mb + ms) 

where Mb^ and fs^ are the mass and decay constant 
of the Bg meson and Nc is the number of colors. B 
and ^5 are defined such that B = Bs = 1 corresponds 
to the factorization (or 'vacuum insertion') approach, 
which can provide a first estimate. Their numerical 
values are known from Lattice QCD calculations. 

The l/rrifc corrections are obtained expanding am- 
plitude Eq. in terms of light quark momentum 
and matching it to four-quark operators that contain 
derivatives [J Q, 

The Si/jn term can be written in the following form: 

6i/m = {{1 + 2z) [cUR2 + 2R4) - 

^2 



- 2 (2C1C2 
12z2 



NcCi) (i?i 



R2) - 



1 -4z 



[(2C1C2 + N,C^) {R2 + 2i?3) 



(23) 
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where additional operators that contain derivatives 
appear 



Ri = 

R2 = 

i?3 = 

i?4 = 



rrib 



^6,Sp(l - -i^)'^'s, 5, (1 - 75)5, (24) 



Their matrix elements are P, Q : 



Bs\R2\Bs) 
{Bs\R3\B,) 
(Bs\R4\Bs) 



= 2 



1 



= -1 



1 



nib 



27V, ) ■'s, £ 



(25) 
l^l B ' 



-I S. 



fl Ml I 



-I B 



Among these S-parameters, Bf and i3| are the most 
widely studied and well known in lattice and light cone 
QCD. In this paper we use the results of Ref. Q. The 
rest of " bag" parameters is estimated using a vacuum 
insertion approximation. The color-rearranged opera- 
tors Ri were eliminated using Fierz identities and the 
equations of motion as in Eq. 

As it was mentioned earlier 0(1/to;,) corrections are 
quite large [l|, Computing 0(l/m^) we directly 
control convergence of 1/m expansion in lifetime dif- 
ference calculation. At this order we get more oper- 
ators that contribute to the Ar{Bs). There are two 
different types of operators. One class of them in- 
volves operators computed by further expansion of 
Eq. ^ - they are called kinetic corrections. An- 
other type arises from interaction of quarks with back- 
ground gluon field. 

The kinetic corrections can be written as: 



24^2 



r(3-10z) [CfWs- 



'^^/"^ - (l-4z)2' 
{2CiC2+N,Ci)iW3 + W2/2)] 

— — 1 [ClQs - (2 C1C2 + N,Cl){Qs + Q/2)] 

1 — 4z mf 



1 -4z 



Wi 
VF4 



rrib 

^5,S"SV(1 - 75 



-\bi'T)°"T)'^ {I - ^50oJ^pSi bj{l - 75)sj , 

- 75)«^M^cS^ hl^i^ - 75)sj , (27) 



where, as before, we have eliminated the color- 
rearranged operators Wi in favor of the operators Wi. 
Due to absence of results from lattice and light cone 
QCD, the parametrization of the matrix elements of 
these operators is given. In the pure factorization ap- 
proach all the bag parameters a; should be set to 1: 



{Bs\Qh\Bs 
(^s\Wx\B, 

(^s\W3\B,) - 

(S,|W4|Bs) = 



2 2 



{nib + 



(28) 



= 1 



1 

2Nr 



Ml 



1^2, 



-1 



1 



r2 X 2 

., ^ / ivi 



Ibs^-'b, 



Ml 



2 JbJ^^b^ 



ml 



Ml 



as 



1 "3 , 



1 a4 , 



In addition to the set of kinetic corrections considered 




[2CfWi - 2 (2 C1C2 + NrCi)iWi + W2/2) 



(1 - 2z)^{Cf + 2C1C2 + N,Cl)Qi 



(26) 



The operators in Eq. ([26)1 are defined as 
Qr = ibiSi)s+p{bj'^3)s+P , 



FIG. 1: Diagrams contributing to corrections due to inter- 
action with background gluon field. 

above, the effects of the interactions of the interme- 
diate quarks with background gluon fields should also 
be included at this order. The contribution of those 
operators can be computed from the diagram on Fig[l] 
resulting in 
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- (1 - 4z) + 4zP3 - 4zP4] 

The local four-quark operators contributing to this 
correction are given in Eq. (|29p . 



Pi 

Pa 

Pa 
P5 
Pe 
Pi 
Ps 



- 75)s, hY{l - l^KiGl^-^i (29) 
-75)5, fefe7.(l-75)ifciG;i.Si, 



bkD D 7^(l-75)t^iG^.s, 6,7o(l-75)s., 

75)Sfc tfeiG;^^ &,7a(l-75)s,, 



1 - i — i^i — a 

—hD D 7'^ 1 



1 - ^ — U < — CK ^ , 



75)sfc ifciG" 6,7,.(l-75)s, 



Following [6(] these operators are parametrized the fol- 
lowing way: 



1 „ „ [Ml \^ 
{B,\P^\Bs) = jflMl ~ Ij A 



(30) 



It is hard to obtain precise prediction for lifctimc dif- 
ference with so many operators contributing. Never- 
theless contribution from 61 /m and 61 can be eval- 



uated. In our numerical calculations we assume the 
pole mass of b-quark to be rrib = 4.8 ± 0.2GeV and 
fs = 230 ± 25MeV. In order to see the effect of 
0(l/m^) corrections we fix all perturbative param- 
eters in the middle of their allowed ranges to show 
dependence of AF^^ on non perturbative parameters 
B,, aj, /3., defined in EqEl [SSI [Ml EQ] 



AFb^ = [0.0005B + 0.1732B, + 0.0024Bi 
-0.0237^2 - 0.0024^3 - 0.0436^4 
+2 X 10"^ai 4 X 10"^a2 + 4 x lO^'^aa 
-h0.0009a4 - 0.0007a5 
-h0.0002/3i - 0.0002/32 + 6 x 10'^ P3 
-6 X IQ-^Pi - 1 X 10-^/35 
-1 X IQ-^Pe + 1 X 10-^/37 -I- 1 X 10-^/?8] (ps"^) 



(31) 



It is obvious that 0(l/m^) corrections provide minor 
effect on calculation of Bg — Bg lifetime difference. 
Contribution from interaction with background giuon 
field is essentially negligible. 

To obtain full SM prediction of AF^^ we vary val- 
ues of parameters of matrix elements. We generate 
100000-point probability distribution of the lifetime 
difference obtained by randomly varying our param- 
eters within ±30% range around their factorization 



value or within zblcr for parameters known from ex- 
perimental data or lattice QCD calculations. The re- 
sulting distribution is presented on the Figl2] There is 




0.03 0.06 0.09 0.12 0.15 0.18 AF, ps" 



FIG. 2: Histogram showing the random distribution 
around the central values of various parameters contribut- 
ing to Bs -lifetime difference AFs^. 

no theoretically consistent way to treat this diagram 
since it is not expected for theoretical predictions to 
have Gaussian distribution. Nevertheless we can give 
a numerical prediction estimating position of peak as 
the most probable value and the peak width at half of 
height as theoretical uncertainty. 



AFs, 
AFs^ 



0.072 ±0-«^4^^- 



0.030 
0.104 ±0.049 



ps 



(32) 



where in latter result we added theoretical error for 
our calculation of AF and experimental error from 
determination of F^^ in quadrature. Additional im- 
provement in lattice or QCD sum rules determination 
of "bag" parameters would make this prediction even 
more solid. 



III. NEW PHYSICS CONTRIBUTIONS TO 
LIFETIME DIFFERENCE 



In the previous section it was shown that 
0(l/TO^)corrections to the lifetime difference of Bs 
and Bs mesons are small. Since we have reliable pre- 
diction of AFb^ it might be interesting to consider 
possible effects of physics beyond the Standard Model 
on the lifetime dilference in B^ system. 

As was pointed out long time ago 0, [1], CP- 
violating contributions to M12 must reduce the life- 
time diffence in iJ^-system, as 



AF, = AFf^^cos^ 29 



(33) 
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where 9s is a CP-violating phase of M12, which is 
thought to be dominated by some A_B — 2 New 
Physics. On other hand, CP conserving AS = 1 
NP amphtude can interfere with SM contribution 
constructively or destructively, depending on the NP 
model. 

There was no spectacular NP phases observed in 
Bs mixing, thus it is important to estimate the CP- 
conserving contribution to AFj. We shall consider it 
using the generic set of effective operators, and then 
apply our results to popular extensions of the SM. 
Using the completeness relation the NP contribu- 



tion to the lifetime difference becomes 



y 
r 



(SJlmTlS, 



^\d''xT{7^f,r\x)7^f^'i^)) 



(34) 



We represent the generic NP A6 = 1 hamiltonian as 
Kp^'^ = E Pi(m)Qi +C2{p)Q2] (35) 



9:9 



where ri_2 are arbitrary combinations of Dirac ma- 
trices and Ci,2(/u) are Wilson coefficients evaluated at 
energy scale /i. This gives us the following contribu- 
tion to lifetime difference: 



Mf 



(36) 



99' 



where z, (3, g, 7 stand for color indices, operators Oj^''^ 
are the following: 

of"'' = {kT'^pr2s^) {hTipr^sp) 

of"" = (foJ-Pas^) (fofeFipr^s^) (37) 
of"" = {kr''pT2S^) {hV^T^sp) 

of"-' = (6,r''r2s^) (6fcrir,s^) , 

with p being a 6-quark momentum, and Ki are the 
following combinations of Wilson coefficients 

Ki - (C2C2iV, + (C2C1 + C1C2)) 

K2 = CiCi (38) 

with the number of colors Nc—3. 

Defining z = ra^/ral and z' = m'^ /ml coefficients 
Ij{z, z') can be written as follows: 



l2{z,z') 

h(z,z') 
h{z,z') 



24mr7r 



[l + {z + z')~2{z-z'f 



^^[l + z'-z] 
Svr 

$ 

47r 



(39) 



z' \1- z' + z] 



where $ is available phase space of process $ = 
m,/2(l-2(z + z') + (z-z')2)'^' 

A. Multi-Higgs model 

One of possible realizations of New Physics is a 
charged Higgs doublet model proposed in 10]. This 
model provides new flavor changing interaction medi- 
ated by charged Higgs bosons, ft leads to the following 
four-fermion interaction: 



AB = 1 
ChH 



b.Fig, q^T2Sj 



(40) 



H 



where P; 



1, 2 are 



Pi = mbV*^cot f3PL — mcV*i, tan 

^2 = TTisVcs cot PPr - nicVcs tan PPl 



(41) 



This leads to three operators with various coefficients, 
matrix elements of which contribute to the ychH- 

ychH 



MbVb aTJ 

[(Qi) {AK2Xsh coe (3 + 2(cot2 Pmlxsh - mbxh){K2 - Ki)) + 
+ (Q2> {-2KiXsh cot2 (3 + (cot^ Pmlxsh - mbxh){K2 - Ki)) + 



+ (Q3)(ifi + K2) (x2 tan2 (3h - m^xh)] 



(42) 



li and Ki were defined above, x — mdmi, and Xs 
nis/mi,. (Qi) are as follows: 



{nib + 



Ql = {bi)L{si)B{bk)R{sk)L 

(Qi) 

Q2 
{Q2) 
Q3 
(Q3) 



2 + — 



{bi)RY{si)ri{bk)Llu{sk)i 



2'" " \ Nc 

{bi)Ll''{Si)Libk)Ll^{sk)L 



(43) 



(44) 



(45) 



1 



flMl { 1 



1 



0.05 



For values of AIh — 856*6^ and cot /3 
gives UchH ~ 0.0034. This is about 10% of Standard 
Model value. Dependence of ychH on mass of Higgs 
boson is given on Figl3] 
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100 105 110 M (GeV) 




FIG. 3: Dependence of ychH on mass of Higgs boson: 
solid line - tan /3 = 20, dashed line - tan P = 10, dotted line 
- tan /9 = 5, dash-dotted line - tan /3 = 3 



(a)Solid lineinon manifest LRSM (k = 1), dashed line; 
manifest LRSM 



B. Left-Right Models 

One of the possible extensions of the SM is a Left- 
Right Symmetric Model (LRSM) which assumes the 
extended SU{2)l x SU[2)h symmetry of the theory. 
In this model additional flavor changing interaction is 
provided by mediating right-handed M^'^'^^ -bosons. In 
this case flavor mixing is described by right-handed 
CKM matrix K-i^^ and 



ri^2 = YPb^ 



^qq — ^cb ^cs 



V2 



(46) 
(47) 



here — g|j./8M^(^^ and for future calculations 
we take gL = kgn 

Such model gives us the following prediction for 
value of y: 




1000 1100 1200 1300 M,., 



(b)Solid line: k = 1, dashed line: k = 1.5, dotted line: k = 2 



FIG. 4: Contribution to ATbs/^Bs in the left-right sym- 
metric models, (a) - Dependence on (b) - Depen- 
dence on Af^' in non-manifest LRSM. 



Ulr 



X Ci(Q2) +^2(02) 



(48) 



One of possible realizations of such scenario which 
gives the biggest numerical value of yLR is a "Non- 
f,^^ w 1) with My 

of Vlr 

LR" {{V}f^ = V^j)) contribut ion from this model is 

less. Dependence of yLR. on mass of VK^^' boson is 
given on Fig. [4] 



manifest LR" (V^/"' w 1) with M^iR) = 1 TeV value 
-0.015 was obtained. In case of "manifest 



IV. CONCLUSIONS 



We computed the subleading 0(l/m^) corrections 
to the lifetime difference of mesons. The correc- 
tions depend on 13 non-perturbative parameters ai 
and Pi. We generated probability distribution of life- 
time difference by varying parameters ±30% around 
their "factorization" values or within la for param- 
eters know from Lattice QCD. The results are pre- 
sented on Figl2l Translating this diagram into numer- 
ical prediction for AFb^/Fb^ we obtained the most 
precise available today theoretical prediction for life- 
time difference: 
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AFb^ 0.072 ± 



0.034 
0.030 



ps 



Ar, 



= 0.104 ±0.049 



(49) 



The effect of 1/ corrections to the hfctime differ- 
ence is small. 

The generic AS = 1 New Physics contribution to 
the lifetime diffcrncc in Bs system is considered. We 



considered four-fermion effective Hamiltionan of the 
generic Standard Model extension and computed its 
contribution to the AF ■ It can reduce or increase 
the SM contribution depending or particular choice of 
the model. Two models of physics beyond the Stan- 
dard Model considered. Contribution of charged Hig- 
gses to the lifetime difference is negligible. LRSM con- 
tribution is significant and parameters of this model 
can be constrained based on AF^^ measurements. 
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